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The dimensionless ratio/ = XM /rjC w relating the thermal conductivity, molecular weight, viscosity, 
and constant volume molar heat capacity has been determined for several nonpolar polyatomic gases 
in the neighborhood of room temperature (270°-295° K). The experimental method, due to Eckert 
and Irvine, provides a direct determination of / by measurement of the subsonic temperature recovery 
factor. A recent theory of Mason and Monchick has been used to calculate collision numbers for 
rotational relaxation from the experimental data as follows: CH 4 , 9.4; CF 4 , 3.0; SF 6 , 2.5; C 2 FL. 2.4; 
C 2 H 6 , 4.0; O 2 , 12; N* f 7.3; CCb, 2.4; and C 2 H 2 , 1.8. Collision numbers for the near-spherical molecules 
were in close accord with a classical theory for rough sphere molecules with attractive forces; ethyl- 
ene, which deviates appreciably from spherical symmetry, exhibited a smaller collision number. The 
data on linear molecules were in qualitative agreement with a quantum treatment. In general, col- 
lision numbers for rotational relaxation arc determined bv the following factors: (1) The molecular 
mass distribution, (2) the strength of the intermolecular attractive forces, and (3) the molecular 

Dor me. 


INTRODUCTION 

I N an earlier paper 1 we reported experimental 
values of the dimensionless ratio 

/ = XAf/,C„ (1) 

which relates the thermal conductivity, molecular 
weight, viscosity, and constant volume molar heat 
capacity. The experimental method, due to Eckert 
and Irvine, 2 provides a direct determination of / 
by measurement of the subsonic temperature re- 
covery factor. Helium, argon, nitrogen, oxygen, and 
hydrogen were studied in the range 90°-290° K. 

The data on helium and argon were in close ac- 
cord with the rigorous Chapman-Enskog kinetic 
theory for monatomic gases, 3 for which / is very 
slightly larger than 2.5. The results for nitrogen, 

1 C. O’Neal, Jr., and R. S. Brokaw, PhvsT Fluids 5, 567 
(1962). 

2 E. R. G. Eckert and T. F. Irvine, Jr., J. Appl. Mech. 24, 
25 (1957). 

3 S. Chapman and T. G. Cowling, The Mathematical 
Theory of Non-Uniform Gases (Cambridge University Press, 
New York, 1952), p. 235. 


oxygen, and hydrogen, as well as Novotny and 
Irvine’s data on carbon dioxide, 4 were analyzed in 
terms of Mason and Monchick’s 5 theory for poly- 
atomic gases, which relates deviations from the 
modified Eucken approximation to the relaxation 
times of internal energy states. Collision numbers 
for rotational relaxation for nitrogen, oxygen, and 
carbon dioxide were found to be nearly independent 
of temperature. The results confirmed the fact that 
the interchange of translational and rotational energy 
is unusually difficult for hydrogen. 

Mason and Monchick’s 5 theory for nonpolar poly- 
atomic gases may be written 



4 L. Novotny and T. F. Irvine, Jr., J. Heat Transfer 83, 
125 (1961). 

5 E. A. Mason and L. Monchick, J. Chem. Phys. 36, 1622 
(1962). 
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Here p is the density and D is the average coeffi- 
cient for the diffusion of internal energy; for nonpolar 
gases it can be taken to be the self-diffusion coeffi- 
cient. Cp trans and C v in fc are the translational and 
internal contributions to the heat capacity C v ; r k 
is the relaxation time for the fcth internal energy 
mode, and C vk is the heat capacity of that mode, 
while p is the pressure. 

The first two terms in Eq. (2) are simply the 
modified Eueken approximation, developed by Chap- 
man and Cowling, 6 Schafer, 7 and Hirschfelder. 8 The 
third term is important only for small relaxation 
times. In small, rigid polyatomic molecules these 
are the rotational relaxation times; in flexible 
molecules, the vibrational relaxation times may also 
be small. (The only non rigid molecule considered 
in this paper is ethane.) 

It is further convenient to express the relaxation 
times in terms of a collision number 


a _ 4 pr k 

— — 

'Too II 7T 7} 


(3) 


Consequently, for nonpolar gases 

f = f m 

/ /mod. Eueken / y \ O / * J *7 V 1 / 

7r(-,, \Z 7} / k /J k 

This paper reports / values in the vicinity of 
room temperature for several nonpolar polyatomic 
molecules, again determined using Eckert and 
Irvine’s 2 experimental technique. Nonlinear mole- 
cules, possessing three rotational degrees of freedom, 
were methane, tetrafluoromethane, sulfur hexafluo- 
ride, ethane, and ethylene. Linear molecules, with 
two rotational degrees of freedom, included oxygen, 
nitrogen, carbon dioxide, and acetylene (some pre- 
liminary results on methane and the linear mole- 
cules have already been reported elsewhere 9 10 11 ). Col- 
lision numbers for rotational relaxation have been 
computed from Eq. (4) and interpreted in terms of 
the existing theories of rotational relaxation. 10,11 
This provides considerable insight as to the molecular 
parameters which determine rotational relaxation 
rates. 


EXPERIMENTAL 


tors, 1 and the method has been discussed at length 
by Irvine. 12 Consequently, the technique will merely 
be outlined here. The method is based on the rela- 
tion describing the adiabatic temperature attained 
by a flat plate in a subsonic gas stream. This re- 
covery temperature T r is conveniently described by 
a recovery factor 

r = (T r ~ r.)/( T t - T m ) t (5) 

where T a is the static temperature of the stream at a 
sufficient distance from the plate and T t is the total 
temperature (the temperature measured at a stagna- 
tion point). Pohlhausen, 13 by integrating the laminar 
boundary layer equations, found that the recovery 
factor is a function of Prandtl number C p rj/M\ alone. 
Spence 14 has shown that for Prandtl numbers near 
unity 


r ~ (C.Tt/m)* = (Cjfcy. (6) 


Indeed Eq. (6) reproduces numerical calculations to 
better than sixteen parts in ten thousand for Prandtl 
numbers between § and l. 9 

To measure the recovery factor, the gas is ex- 
panded through a convergent nozzle; the pressure 
ratio is such that the issuing gas remains slightly 
subsonic. A butt-welded differential thermocouple is 
suspended along the nozzle axis with one junction 
in the low velocity region upstream of the nozzle 
sensing the total temperature and the other junc- 
tion just downstream of the nozzle exit assuming the 
flat plate recovery temperature. The static tempera- 
ture is not measured; rather, it is computed from the 
pressure ratio across the nozzle assuming iscntropic 
flow. 

The dimensionless entropy of a gas may be 
written 15 


S 

R 


S° i 

R~ lnp - 


p dB 
RdT + 


(7) 


Here S is the entropy, S° is the entropy of the ideal 
gas at one atmosphere, R is the gas constant, p 
is the pressure, B is the second virial coefficient , and 
T is the absolute temperature. Consequently, for 
isentropic flow, 


We have previously described the apparatus and 
experimental procedure for measuring recovery fac- 



C[0 

R 


-In 


Pt _ Pt ( dB\ 
p s R \dT) 


T -T t 




• ( 8 ) 


6 Reference 3, p. 238. 

7 K. Schafer, Z. Physik. Chem. B53, 149 (1943). 

*J. O. Hirschfelder, J. Chem. Phys. 26, 282 (1957) 

9 R. S. Brokaw and C. O’Neal, Jr., Ninth Symposium 
{International) on Combustion (Academic Press Inc., New 
York, 1963), p. 725. 

10 N. F. Sather and J. S. Dahler, J. Chem. Phvs. 37, 1947 
(1962). 

11 R, Brout, J. Chem. Phys. 22, 1189 (1954). 


The last two terms in Eq. (8) are a small correction 


12 T. F. Irvine, Jr., Ph.D. Thesis, University of Minnesota, 
1956. 

13 E. Pohlhausen, Z. Angew. Math. u. Mech. 1, 115 (1921). 

14 D. A. Spence, J. Aerospace Sci. 27, 878 (1960). 

15 J. O. Hirschfelder, C. F. Curtiss, and R. B. Bird, 
Molecular Theory of Gases and Liquids (John Wiley & Sons, 
Inc., New r York, 1954j_ Ec[. 3. B-3, p. 231, 




THERMAL CONDUCTIVITY AND VISCOSITY FOR NONPOLAR GASES 1677 


due to deviations from the ideal gas law. To com- 
pute the stream temperature an approximate value 
was obtained from Eq. (8) using tabulated standard 
entropies and the experimental pressure ratio. The 
derivatives of the virial coefficients were then cal- 
culated assuming a Lennard-Jones (12-6) potential 
together with appropriate force constants. 16 Equa- 
tion (8) was then used to calculate the final value of 
S° t ./R , and the accurate stream temperature was 
obtained by interpolation in the table of entropies. 
Recovery factors and values of / could then be calcu- 
lated from Eqs. (5) and (6). Entropies and heat 
capacities for oxygen, nitrogen, and ethane were 
taken from National Bureau of Standards compila- 
tions, 17 ’ 18 while values for the remaining species were 
calculated by the methods discussed by McBride, 
Heimel, Ehlers, and Gordon, 19 using their spectro- 
scopic constants and computer program. 

The present measurements were all carried out at 
room temperature; hence, temperature control was 
no problem. Furthermore, it was possible to work 
with fairly high flow rates so that the pressure ratios 
and differential thermocouple emfs could be read 
with greater accuracy. As a consequence, the present 
measurements show considerably less scatter than 
our previous data. 1 

In order to detect instrumental errors, 83 check de- 
terminations of * the recovery factor of argon were 
interspersed throughout the course of this work. 
These values correspond to a temperature of about 
275° K and yield a mean value of / = 2.5028 with a 
standard error of 0.0015. (The standard deviation 
of each datum was 0.0135.) This is in excellent agree- 
ment with the theoretical value of f = 2.5023, and 
may be regarded as a calibration of the apparatus. 

RESULTS AND DISCUSSION 

Experimental / values are presented in Table I 
together with the corresponding temperatures. Also 
shown are values of Z~ 0 \, the reciprocal of the col- 
lision number for rotational relaxation, calculated 
from Eq. (4). In calculating Z~\ the relation pD/t\ = 
(f )A* has been invoked. 20 A* is a very slowly varying 
function of temperature; values were obtained assum- 
ing either the exponential —6 or Lennard- Jones 
(12-6) potential. Also shown in Table I are mean 


temperatures T; reciprocal collision numbers Z~ Q \, 
standard deviations a; and standard errors of the 
mean e. 

For all of the molecules except ethane a thousand 
or more collisions are required for vibrational relaxa- 
tion; hence, vibrational effects could be ignored in 
Eq. (4) . In ethane the mode corresponding to 
hindered rotation about the carbon-carbon bond 
has a relaxation time of 1.24 X 10“ 9 sec, or 17 col- 
lisions, at 296.4° K. 21 Assuming that a 3 kcal barrier 
hinders internal rotation, the tabulations of Pitzer 
and Gwinn 22 give C/R = 0.999 for this mode at 
286° K. The rotational collision numbers for ethane 
have been calculated assuming this heat capacity 
and collision number for internal rotation; if the 
internal rotation were ignored, ethane Z~ Q \ values 
would all be larger by 0.040. 

Collision numbers determined from recovery 
factor, or /, are compared with other measurements 
in Table II. The present values are generally in 
accord with those obtained by other techniques, 
within the admittedly rather large uncertainties asso- 
ciated with such determinations. The results for 
oxygen seem to fall into two groups: collision num- 
bers in the ranges of 2-6 and 12-30. Our value falls 
in the latter category. The only other serious dis- 
agreement is in the case of carbon dioxide, where 
acoustical measurements indicate that sixteen col- 
lision s are required for relaxation. We are inclined 
to believe the acoustic result is in error — we can 
think of no property of the carbon dioxide mole- 
cule which would cause such long relaxation times. 

In order to compare our collision numbers with 
theoretical calculations, it is desirable to divide the 
molecules into two groups, linear and nonlinear, 
possessing two and three rotational degrees of free- 
dom, respectively. We chose to discuss the nonlinear 
molecules first, since a classical theory for nonlinear 
molecules with attractive intermolecular forces 
(rough spheres and spherocylinders surrounded by 
square wells) has recently been developed. 10 

Nonlinear Molecules 

Sather and Dahler 10 obtained the following ex- 
pression for the rotational relaxation time of a rough 
sphere with a spherically symmetric, square, at- 
tractive well: 


16 Reference 15, Tables I-A and I-B, pp. 1100-1115. 

17 J. Hilsenrath, et «/., Tables of Thermal Properties of 
Gases, Natl. Bur. Std. Circ. 564, (1955). 

18 F. D. Rossini, et ah, Selected Values of Properties of 
Hydrocarbons , Natl. Bur. Std. Circ*. 461 (1947). 

19 B. J. McBride, S. Heimel, J. G. Ehlers, and S. Gordon, 
NASA Special publication SP-3001, 1963. 

20 Reference 15, Eq. 8.2-48, p. 540. 


Tro t 


_ 16 2 (4//W) 

3 m [1 + (4 1/ma 2 )} 2 



(9) 


21 L. M. Valley and S. Legvold, J. Chem. Phys. 33, 627 
(1960). 

22 Iv. S. Pitzer and W. D. Gwinn, J, Chem. Phvs. 10, 
428 (1942). 
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Table I. Experimental Results 


T° K / Z" 1 T°K f Z~ l T°K f Z 1 T° K / Z -i 


Methane 


278.5 

1.813 

0.141 

281.9 

1.807 

0.152 

279.0 

1.831 

0.098 

282.0 

1.846 

0.056 

280.0 

1.844 

0.064 

282.3 

1.812 

0.137 

280.0 

1.819 

0.125 

282.8 

1.811 

0.140 

280.0 

1.814 

0.136 

283.0 

1.854 

0.037 

280.2 

1.821 

0.118 

283.3 

1.810 

0.135 

280.8 

1.821 

0.118 

284.0 

1.857 

0.028 

281.0 

1.834 

0.087 

284.1 

1.824 

0.106 

281.3 

1.822 

0.117 

284.1 

1.821 

0.114 

281.5 

1.822 

0.117 

285.4 

1.828 

0.095 

T = 

282 °K; 2r 

hot = 0. 

.106; 9 = 

0.039; e 

= 0.009 



Carbon tetrafluoride 


286.9 

1.527 

0.328 

289.1 

1.522 

0.351 

287.3 

1.528 

0.325 

289.6 

1.529 

0.313 

287.7 

1.526 

0.334 




T = 

288 °K; Zr 

hot = 0, 

330; 9 = 

0.014; e 

= 0.006 



Sulfur Hexafluoride 


288.3 

1.454 

0.284 

291.2 

1.444 

0.354 

289.2 

1.440 

0.391 

291.4 

1.437 

0.407 

290.5 

1.432 

0.445 

291.5 

1.445 

0.343 

291.0 

1.453 

0.283 

291.6 

1.437 

0.403 

291.1 

1.449 

0.316 




T = 

291 °K; Zr 

hot - 0. 

358; 9 = 

0.057; e 

= 0.019 



Ethylene 



281.2 

1.624 

0.410 

282.6 

1.622 

0.412 

281.9 

1.616 

0.434 

282.8 

1.621 

0.417 

282.2 

1.620 

0.421 

286.6 

1.613 

0.430 

282.2 

1.617 

0.430 

287.7 

1.615 

0.420 

T = 

283 Tv; Zr 

hot = 0. 

422; 9 = 

0.009; e 

= 0.003 



Ethane 



284.2 

1 .585 

0.265 

286.3 

1.578 

0.286 

284.9 

1.596 

0.223 

287.3 

1.593 

0.224 

285.4 

1 .578 

0.288 

288.0 

1.594 

0.217 

285.7 

1.585 

0.262 

288.0 

1.594 

0.217 

T = 

286 Tv; 2r 

'h ot - 0.248; 9 = 

0.031; e 

= 0.011 



Oxygen 



276.1 

1.988 

0.076 

280.7 

1.967 

0.134 

276.6 

1.991 

0.069 

280.8 

1.989 

0.074 

277.3 

2.001 

0.041 

281.0 

2.002 

0.036 

277.3 

1.971 

0.124 

281.3 

1.962 

0.148 

278.2 

2.003 

0.036 

281.6 

1.990 

0.069 


Oxygen 


278.3 

1.978 

0.105 

282.0 

1.961 

0.150 

278.3 

1.982 

0.093 

282.0 

1.987 

0.078 

279.4 

1.991 

0.069 

282.4 

1.987 

0.079 

279.5 

1.995 

0.057 

282.6 

1.967 

0.134 

279.9 

2.000 

0.043 

282.8 

1.966 

0.136 

280.4 

1.990 

0.069 

282.9 

1.991 

0.067 

f = 

280 °K; Z~ ] 

hot - 0. 

086; a = 

0.037; e 

= 0.008 



Nitrogen 



274.3 

1.977 

0.122 

280.7 

1.984 

0.102 

274.3 

1.963 

0.160 

281.1 

1.974 

0.130 

275.4 

1.974 

0.129 

281.5 

1.959 

0.171 

275.6 

1.974 

0.129 

281.7 

1.972 

0.134 

275.9 

1.979 

0.117 

281.8 

1.961 

0.165 

276.4 

1.978 

0.118 

282.3 

1.969 

0.144 

277.0 

1.986 

0.097 

282.4 

1.973 

0.133 

277.3 

1.962 

0. 162 

282.6 

1.966 

0.152 

277.5 

1.981 

0.109 

282.7 

1.967 

0.148 

278.8 

1.988 

0.092 

282.9 

1.968 

0.147 

278.9 

1.978 

0.119 

283.3 

1.961 

0.165 

279.6 

1.960 

0.169 




T = 

279 °K; Z~ 

hot = 0. 

138; 9 = 

0.022; e 

= 0.005 



Carbon Dioxide 



282.0 

1.724 

0.454 

284.0 

1.740 

0.388 

282.0 

1.723 

0.457 

286.4 

1.729 

0.425 

283.0 

1 .732 

0.420 

287.4 

1.731 

0.410 

283.0 

1.731 

0.425 

288.1 

1.731 

0.408 

284.0 

1 . 732 

0.417 




T = 

284 °K; Z 

hot - 0. 

II 

lb 

W 

0.021; e 

= 0.007 



Acetylene 3 



283.3 

1.649 

0.457 

288.6 

1.624 

0.558 

283.6 

1.662 

0.394 

288.6 

1.612 

0.614 

284.0 

1.662 

0.395 

288.8 

1.602 

0.663 

284.9 

1.660 

0.399 

289.1 

1.651 

0.426 

285.0 

1.661 

0.395 

289.4 

1.613 

0.608 

286.0 

1.609 

0.638 




286.5 

1.606 

0.650 

289.6 

1.608 

0.629 

286.9 

1.609 

0.635 

290.0 

1.607 

0.633 

287.3 

1.610 

0.629 

290.2 

1.610 

0.619 

287.8 

1.608 

0.638 

290.7 

1.656 

0.395 

T = 

287 Tv; Zr 

hot = 0, 

,546; 9 = 

0.111; e 

= 0.025 


» The data on acetylene fall into two distinct groups. For the most recent data Z~ l ro t “ 0.626, whereas, for the earlier results Z~ l ro t = 0.409. 
This is the cause of the high standard deviation We have no valid reason to reject the earlier data although we have been unable to reproduce them. 


Here n is the number of molecules per cm 3 , I is the 
moment of inertia, m is the molecular mass, k is 
the Boltzmann constant, a is the diameter of the 
rough sphere core, while g(cr) is the value of the radial 
distribution function at a. In the low density limit, 23 
g{a) — exp (e/fcT 7 ), where e is the depth of the well 
at cr. 

Two features of Eq. (9) are noteworthy. First, 
since the quantity 47/mo- 2 is generally less than 0.2, 


23 Reference 15, Eq. 4.9-1, p. 321, 


the relaxation time is only weakly dependent on a , 
the position of repulsive core. Secondly, since there 
is no rotational energy transfer accompanying the 
velocity impulse at the outer edge of the potential 
well, the only contribution to the rotational relaxa- 
tion stems from the impulse at the core. Thus, the 
width of the well is not important; in fact, Eq. (9) 
should apply to a Sutherland potential with a rough 
core — that is, a potential with an inverse sixth power 
attractive portion to account for van der Waals 
forces. 
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Table II. Comparison of Rotational Collision Numbers 


Gas 


Experimental method 




Recovery 

a Acous- b Low pressure 

Shock 

Impact 


factor 

tical thermal 

conductivity 

thick- 

ness 

tube 

n 2 

7 

5.3, 6, 4-6 

5.5 C 

<7, d 

<14® 

Oo 

12 

2-4,4.1,12, 20 

14,12-30 

7 a 


h 2 

Large** 

240-360 300 

>150 f 

160, d 
310® 

co 2 

CH< 

2.4 

9 

16 

14-17 

1-2® 



• Cited by T. L. Cottrell and J, C. McCoubrey, Molecular Energy 
Transfer in Gases (Butterworths Scientific Publications Ltd., London 
1961) Ch. f>. 

b F. G. Waelbroeck and P. Zukerbrodt, J. Chem. Phys. 28, 524 (1958). 
' W. H. Andersen and D. F. Hornig. Mol. Phys. 2, 49 (1959). 
d P. W. Huber and A. Kantrowitz, J. Chem. Phys. 15, 275 (1947). 
« W. Griffith, J. Appl. Phys. 21, 1319 (1950). 

! E. F. Greene and D. F. Hornig. J. Chem. Phys. 21, G17 (1953). 

« Reference 1. 


The coefficient of viscosity is 

V = (5/16o-^)(mkT/T) i , 


where is the square of the viscosity collision 
diameter [ = <r 2 Q (2,2> * in the notation of reference 15, 
page 527, Eq. 8.2-10]. Thus, from Eqs. (3) and (9) 
together with the viscosity formula and the perfect 
gas law, 


Z~\ = 


r (4 1/mo-l) l _e_\ 

12 [1 + (4//W)] 2 \kT) 


f 1 f\ n \ 



(10b) 


We might expect that the rough sphere model with 
an attractive potential well would apply to non- 
linear molecules which are approximately spherical. 

The experimental results on nonlinear molecules 
are analyzed in Table III. Viscosity collision diam- 



eters were calculated assuming a Lemiard-Jones 
(12-6) potential for all molecules except methane; 
for methane the exponential — 6 potential was used. 
Methane, carbon tetrafluoride, and sulfur hexafluo- 
ride are approximately spherical as shown in Fig. 1. 
Indeed, Z~\ values for these molecules calculated 
from both Eqs. (10a) and (10b) are in satisfactory 
agreement with experiment. [The approximation 
(10b) is actually in closer accord with experiment, 
but it is not clear that this is significant.] The mass 
distribution parameter 4 1/ma 2 v varies more than 
threefold between CH 4 and CF 4 ; thus, the large col- 
lision number for methane is undoubtedly a direct 
consequence of the molecule’s small moment of 
inertia. Note, too, that the effect of attractive forces 
[exp ( e/kT )] is appreciable and roughly doubles the 
transition probabilities. 

The Z~\ values for ethylene and ethane were cal- 
culated using average moments of inertia. Values 
for ethane are also in close accord wdth experiment, 
but in the case of ethylene the agreement is not so 
good; the predicted Z~\ are roughly half of experi- 


Table III. Comparison of Theory and Experiment for Nonlinear Molecules 



Methane 

Carbon 

tetrafluoride 

Sulfur 

hexafluoride 

Ethylene 

Ethane 

Type of potential 

exp-6 a 

L-J 12-6 

L— J 12—6 

L-J 12-6 

L— J 12—6 

U/k, °K . 

152.8 

134.0 

222.1 

224.7 

215.7 

Force constants < a or r m , A 

4.206 

4.662 

5.128 

4.163 

4.443 

l « 

14 





Moments of inertia, 

5.341 

147.3 

307.1 

5.75 

11.03 

gm cm 2 X 10 40 




28.09 

42 .28( 2 ) 





33.84 


4I/m<r% 

0.0476 

0.1612 

0.1375 

0.0784 

0.0929 

exp (e/kT) 

1.72 

1.59 

2.14 

2.21 

2.13 

Z~' rot [ Eq. (10a)] 

0.098 

0.249 

0.297 

0.195 

0.217 

Z~\ ot[Eq. (10b)] 

0.107 

0.335 

0.384 

0.227 

0.259 

Z Lot Experiment 

0.106 

0.330 

0.358 

0.422 

0.248 

Standard error of the mean 

0.009 

0.006 

0.019 

0.003 

0.011 


Reference 28. 
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Fig. 2. Experimental and theoretical values of / for carbon 
dioxide. 

ment. From Fig. 1 it is apparent that the ethane 
molecule is approximately spherical, whereas the 
ethylene structure is definitely less compact and 
symmetric. 

It appears, then, that the classical kinetic theory 
for a rough spherical molecule with an attractive 
well [Eq. (10a)] may provide a method for calculat- 
ing a lower limit to the collision probability for ro- 
tational relaxation, Z T ~\, for nonlinear molecules. 
Molecules such as CH 4 , CF 4 , SF 6 , and C 2 H Q are 
reasonably represented by this model. Less sym- 
metric molecules such as C 2 H 4 exhibit larger transi- 
tion probabilities, or shorter relaxation times, per- 
haps due to longer range asymmetries in the inter- 
molecular force field. 

Linear Molecules 

Brout 11 has developed a quantum mechanical 
treatment for the rotational relaxation of homo- 
nuclear diatomic molecules. He idealizes one of the 
molecules as a mass point and expresses the inter- 
molecular potential 4> in the form 

#(r) = Ae~ pr [ 1 + tP 2 ( cos 0)]. (11) 

Here P 2 (cos 6) = |(3 cos 2 6 — 1) is the Legendre 
polynomial of order 2, and 0 is the angle between 
the line of centers of the colliding molecules and the 
figure axis of one of the molecules. \p is a parameter 
describing the deviation of the potential from spheri- 
cal symmetry. For such a potential, the transition 
probability may be written 

Z~\ = Ut 2 (4I/m<T 2 v ). (12) 

Although Brout developed Eq. (12) for homonuclear 
diatomic molecules, it should apply as well to sym- 
metrical linear polyatomic molecules such as carbon 
dioxide, acetylene, etc. Note that the mass distri- 


bution parameter 4 1 /m<r\ is again relevant, as in the 
classical rough sphere case. 

The potential function given in Eq. (11) is entirely 
repulsive, whereas real molecules have long range 
attractive forces. It would seem reasonable to ac- 
count for the attractive forces by a multiplicative 
factor exp ( e/JcT ) by analogy with Eq. (10): 

z ~ - ! *’(S) esp fe)- <l3) 

Equation (13) indicates some variation of collision 
number with temperature. Assuming that \p is inde- 
pendent of temperature — and, as we shall see, 
Brout’s treatment indicates that it is — this variation 
can be calculated. Thus, over the range of experi- 
mental measurement — roughly 100° to 300° K l — 
Eq. (13) indicates that Z~\ for nitrogen should de- 
crease by 25% and Z~ 0 \ for oxygen should decrease 
by 35%. This temperature variation represents the 
somewhat scattered data 1 on these molecules about 
as well as a temperature independent collision num- 
ber. The experimental data on carbon dioxide 4 ex- 
tend from about 280° to 460° K and are more pre- 
cise. Over this range Eq. (13) indicates a 9% de- 
crease in Z~\. In Fig. 2, the carbon dioxide data are 
shown together with calculated curves assuming: 
First, a temperature independent collision number, 
next, the temperature variation of Eq. (13), and 
finally, the temperature variation of Parker’s 24 classi- 
cal two-dimensional treatment for diatomic mole- 
cules. The data are best represented by a constant 
collision number, although Eq. (13) is also reasonably 
satisfactory. Parker’s temperature dependence would 
appear to be ruled out. It might be noted that Sather 
and Dahler 10 calculated the effect of temperature on 
the relaxation times for spherocylinders and found 
a somewhat less pronounced temperature variation 
than in the rough sphere case. Thus, it would seem 
that if there is a long range asymmetry in the at- 
tractive potential there may be less effect of tem- 
perature on the relaxation time. 

The experimental results on linear molecules are 
analyzed in Table IV. Experimental values of P 
have been computed from Eq. (13); it is seen that 
yp is of the order of unity, varying from 0.7 for 
oxygen to 1.7 for acetylene. Brout 25 has shown how \p 
for diatomic molecules can be calculated assuming 
that the total potential is made up of the sum of 
atom-atom interactions, 

Hr) = x 

24 J. G. Parker, Phys. Fluids 2, 449 (1959). 

25 It. Brout, J. Chem. Phys. 22, 934 (1954). 
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Table IV. Comparison of Theory and Experiment for Linear Molecules 


Carbon 



Oxygen 

Nitrogen 

dioxide 

Acetylene 

Tvpe of potential 

“exp— 6 

6 exp — 6 

C L-J 12-6 

C L-J 12-6 

r«/*, °k . 

132 

101.2 

195.2 

231.8 

Force constants < <r or r m , A 

3.726 

4.011 

3.941 

4.033 

l « 

17.0 

17.0 



Moment of inertia, 

19.3 

14.01 

71.47 

23.65 

gm cm 2 X 10 40 
Data analysis 
41 fm<r\ 

0.111 

0.086 

0.189 

0.094 

exp ( e/kT ) 

1.60 

1.44 

1.99 

2.24 

expt, [Eq. (13)] 

0.74 

1.11 

1.12 

1.71 

\p calc I 2 f 

0.55 

0.48 

1.16 

1.48 

1.10 

1.01 

1.63 

1.80 


a \J. T. Yanderslice, E. A. Mason, and W. G. Maisch, J. Chem Phys. 31, 738 (19">9). 
b *Reference 28. 

OL A. Svehla, NASA Technical Report R-132, 1963. 


where <i> } , is of the form 

= A exp (-jdr,y). 

In this event 

*dUt«ni. ^ m 2 /[12 + hm 2 l (14) 

where < 8 is the index of the repulsive potential and d 
is the intemuelear distance. From Eq. (14) it is clear 
that lies in the range zero to two, which is in agree- 
ment with the experimental ^ listed in Table IV. 

Values of fi can be obtained from two sources. 
For the exponential repulsive potential between light 
atoms and molecules involving hydrogen 26 and 
oxygen and nitrogen 27 fi lies in the range 2.5 to 3.5. 
On the other hand, one can deduce values from 
force constants for the exponential —6 potential 
(fi = a/r m ) ; for light molecules /3 so obtained lie 
in the range 3.5 to 4.5. 27,28 

Accordingly, values of i p have been calculated 
from Eq. (14) assuming fi of 2.5 and 4.5 so as to 
cover the extreme range. For carbon dioxide or 
acetylene, d in Eq. (14) wsls taken as the intra- 
molecular oxygen-oxygen or hydrogen-hydrogen dis- 
tance. In Table IV we see the experimental $ for 
oxygen and acetylene lie between the computed ex- 
tremes while the experimental for nitrogen and 
carbon dioxide lie just outside the calculated range. 
In the case of the acetylene and carbon dioxide 
molecules one might expect the effective centers of 
repulsion to lie inside the hydrogen or oxygen atoms. 

It is not difficult to extend Brout’s treatment to 
molecules with multiple centers of repulsion (say 
one for each atom), and also to apply a correction 
for molecules which deviate more severely from 

26 J. T. Vanderslice, S. Weissman, E. A. Mason, and R. J. 
Fallon, Phys. Fluids 5, 155 (1962). 

27 K. S. Yun and E. A. Mason, Phys. Fluids 5, 380 (1962). 

28 Reference 15, Table 3.7-2, p. 181. 


spherical symmetry. However, because the uncer- 
tainty in the fi values remains, there is no substantial 
improvement over the results shown in Table IV. 

It seems, then, that Brout’s theory, corrected for 
attractive forces, describes the rotational relaxation 
of linear molecules, at least approximately. Accurate 
a priori calculation of relaxation times requires an 
accurate knowledge of the molecular phrenology, 
characterized by fid; such information is not gene- 
rally available. 

In a previous paper 9 we noted an empirical corre- 
lation between the relaxation times of linear mole- 
cules and the molecular quadrupole moments. Ac- 
cordingly, w^e attempted to express the potential in 
the form of Eq. (11) w r ith the angularly dependent 
portion arising from the quadrupole-quadrupole in- 
teraction. First, it w-as necessary to average the po- 
tential over all orientations of one of the mole- 
cules. 29 Next, w r e had to drastically approximate a 
cos 4 8 angular dependence as 

cos 4 8 ~ cos 2 8 cos 2 8 — \ cos 2 8 

to obtain the admittedly very approximate ex- 
pression 

if ~ I Q 4 /r 10 kT$. 

Here Q is the quadrupole moment. If we assume that 
transition occurs w r here is of the order of JcT and r 
is of the order of the viscosity collision diameter, 

* ~ (is) 

Consider carbon dioxide as an example: Bucking- 
ham 30 has recently determined the quadrupole mo- 
ment as —8.2 X 10~ 26 esu. From Eq. (15) we ob- 
tain, for a temperature of 288° K, ^ ~ 0.1. Thus, 

29 Reference 15, pp. 27-28, 985. 

30 A. D. Buckingham, Abstracts of Papers, 144th National 
Meeting, American Chemical Society, Los Angeles, California, 
March 31 - April 5, 1963. 
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we suspect that the quadrupole interaction is not 
important in inducing translation-rotation transi- 
tions; rather,* the correlation previously observed 
arises because the deviations from spherical sym- 
metry in the intramolecular charge distribution 
which cause an asymmetry in the repulsive potential 
also give rise to a quadrupole moment. In other 
words, i p and Q have a common cause, but are not 
directly related. This conclusion must remain tenta- 
tive, however. The derivation of Eq. (15) is very 
crude, and the dependence on distance is very strong. 
(If one assumes the zero energy collision diameter 
in place of the viscosity diameter, ^ ~ 0.4.) 

CONCLUSIONS 

Since the data on nitrogen, oxygen, carbon dioxide, 
methane, and hydrogen 1 yield collision numbers for 
rotational relaxation which are generally in accord 
with values obtained by other experimental methods, 
we conclude that the approximate theory of Mason 
and Monchick 5 is substantially correct. 

Rotational relaxation times for near-spherical 
molecules such as methane, carbon tetrafluoride, 
sulfur hexafluoride, and even ethane are in good 
accord with theory for classical rough sphere mole- 
cules with attractive forces — a reasonable model for 
these substances. Ethylene, which deviates markedly 


from spherical symmetry, has a smaller collision 
number — about half that calculated for the rough 
sphere molecule. 

Collision numbers for linear molecules — nitrogen, 
oxygen, carbon dioxide, and acetylene — are in at 
least qualitative accord with an extension of BroutV 1 
quantum theory for rotational relaxation of homo- 
nuclear diatomic molecules. This theory indicates 
that the deviation of the repulsive force field of a 
molecule from spherical symmetry has an important 
effect. 

In conclusion, it appears that the following factors 
are of profound importance in determining rota- 
tional relaxation times for non-polar molecules: 

1. The mass distribution (characterized by the 
parameter 4 1 /via 2 ). 

2. The strength of the intermolecular attractive 
forces (characterized by e/JcT). 

3. The deviation of the molecular force field from 
spherical symmetry. 
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